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Abstract 

In a recent paper, we included two loop, relativistic one loop and second order 
relativistic tree level corrections, plus leading nonperturbative contributions, to 
obtain a calculation of the lower states in the heavy quarkonium spectrum correct 
up to, and including, 0(a 4 ) and leading A 4 /m 4 terms. The results were obtained 
with, in particular, the value of the two loop static coefficient due to Peter; this 
been recently challenged by Schroder. In our previous paper we used Peter's result; 
in the present one we now give results with Schroder's, as this is likely to be the 
correct one. The variation is slight as the value of b\ is only one among the various 
0(a 4 ) contributions. With Schroder's expression we now have, 

m b = 5 00lle6 4 MeV; m b {rht) = 4 440lg MeV, 

m c = 1 866+i|2 MeV; rh c (rh 2 c ) = 1 531+^ MeV . 

Moreover, 

T(r -> e+e~) = 1.07 ± 0.28 keV (exp. = 1.320 ± 0.04 keV) 
and the hyperfine splitting is predicted to be 

M(T) - Af (77) = 47±l3 M eV . 
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1 Introduction 

In a recent paper' 1 !, to be hereafter denoted by PYI, we evaluated the heavy quarkonium spectrum to 
order . The ingredients for the calculation were the one loop corrections' 2 ! to the static potential, and 
its first and second order contributions to the spectrum; the relativistic, and mixed relativistic-one loop' 3 ' 
corrections; the two loop corrections to the static potential, and the leading nonperturbative corrections' 4 !. 
The two loop correction to the static potential used was that calculated by Peter' 5 !; but the calculation of 
Peter has recently been challenged by Schroder' 6 !. In this last reference, the result of Peter's is checked for 
all pieces except one of the contributions to the C\ coefficient, where an error in Peter's evaluation is pointed 
out. In the present paper we give the results of the calculation using the Schroder results. Note that the 
variation is very small; the reason is that the two loop static potential is only one of the several contributions 
to the 0(a 2 ) calculation. For example, Using the value of the two loop coefficient found by Peter, we found 



m b = 5 015^70° M eV; m. 



1884; 222 MeV 



(P) 



(1.1a) 



to which correspond the MS masses, 



m h {m\) = 4453^32 MeV; m c (m 2 ) = I 547j^ MeV . (P) 



With Schroder's result, one now has 



Tfii, = 5 OOlt™ 4 MeV; 
fh b (fh 2 b ) =4440lg MeV (S) 



(1.1b) 



and, for the c quark, 

m c = 1 866±1§2 McV ; 
m c (m 2 ) = 153lt}?,7 MeV (S) 
For the leptonic decay of the Y, and the hyperfine splitting the results with Schroder, s value of b\ are 



(1.1c) 



r(Y -> e+e") = 1.07 ± 0.28 keV (exp. = 1.320 ± 0.04 keV) 



and 

M{T) - M( V ) = 46.6±^;f McV 
they are almost identical to those obtained before. 
2 The Effective Potential 

We follow the method of effective potentials, and the renormalization scheme of ref. 3. The Hamiltonian 
for quarkonium is 

H = H^+H 1 (2.1a) 



where 



m r 



a s (/i 2 ) ^asip 2 ) \ 1 + ( ai + 



7eA)\ « s (m 2 ) 



and will be solved exactly. Hi is 

Hi = Vtrcc + V} L) + V 2 (L) + + V s . rcl + V spin , 



(2.1b) 



(2.1c) 
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and (we repeat the formulas of PYI for ease of reference) 

V - _1 A 2 + Cf<Xs A 

T/ (i) _ -CF/3 a s (^ 2 ) 2 logr/z 
Kl - 2^ ~ 



= -C7 F 4 L) a^ logr/x ( 2 .ld) 



7r^ r 



y(LL) _ -CF/3oas log T\X 

4-7T 2 r 



s.rel 



2rar 2 



^spin = 3m2 «(« + l)<S(r). 

Here the running coupling constant has to be taken to three loops. a± was calculated in ref. 1 and in 
ref. 3. Both have been checked by independent calculations; the only discrepancy lies in the value of the 
coefficient 1 b\. According to Peter, ^ 

_ [iga + f C(3 )] C A T F n f - [f - 16C(3)] C F T F n f + ^T F n}} ("a) 
~ 24.30, 



while Schroder^ 6 ] gives 

~[W + ¥ C(3)] CXTprif [f - 16C(3)] CyZVn, + f>T 2 n 2 } (2.2b) 
~ 13.2. 



Although the difference between the two lies only in one piece of the coefficient of jgCp, a 47r 2 vs. a 
67r 2 , the Schroder result makes the two loop correction much smaller (and hence the overall calculation more 
believable). We present here results with (2.2b) since they are the ones more likely to be correct (moreover, 
results with (2.2a) can be found in PYI). 

It should be noted that Hi contains a tree level velocity correction, and a velocity- dependent one loop 
piece, Ki.rci- This is because the average velocity in a Coulombic potential is ~ a s , hence a calculation 
correct to order a A s requires tree level 0(v 2 ) and one loop 0(H) contributions. All these terms in Hi may 
be treated as perturbations to first order, except . For this, the second order perturbative contribution 
is also required as this also produces a correction of order a A s . It is to be noted that all the dependence on 
bi is contained in the bi dependence of a s , Eq. (2.1b). 

A last comment concerns the renormalization scheme. We have followed ref. 3 in renormalizing a s in the 
MS scheme; but the mass m that appears in Eqs. (2.1) is the two loop pole mass. That is to say, it is defined 
by the equation, 

S^{ r = m,m) = (2.3) 

where Szij/) ,m) is the quark propagator to two loops. One can relate m to the MS parameter, also to two 
loop accuracy, using the results of refs. 7. 

1 For the values of the constants other than &i entering above formulas, cf. PYI. 
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Nonperturbative corrections are not included in Eqs.(2.1); they will be incorporated later. 



3 Energy Shifts, Order a 4 , A i /m A 

Taking into account the expression for the Hamiltonian, Eq. (2.1), we write 



cw 



E nl = 2m- + ]T 5$ } E nl + 5^ L) E nl + S NP E nl . 



(3.1) 



We define generally the analogue of the Bohr radius, 



a( M 2 ) 



and then (PYI), 



4? E nl 



1 



21 + 1 8n 

(3 C F a 2 (u 2 ) 



C F a s 21 + 1- An 



m 2 n 4 (2l + l)a s ' 



2im 2 a 



(1) _ C F c^al 

ir( L ) nl ~ 9 9 



log ^ + V(n + / + 1) 



log^ + ip(n + I + 1) 



c(l) 



C F (3 2 a* 



+^(n + / + l) 2 +^'(n + / + l) 



+ ' ( ^^ 2) r(n + / + l) g i!(n-Z-i-l) 



l) 2 J' 



c(l) 



V 8 . rel - C 'ni 



E n l = 



m n 3 (2; + l)a 2 ' 

For the vector states (T, Y' , Y"; J/ip, ip' , . . .) one has to add the hyperfine shift, at tree level, 



li) 



3n 3 m 2 a 3 ' 



(3.2a) 
(3.2b) 
(3.2c) 

(3. 2d) 

(3.2e) 
(3.2f) 



The calculation of the second order contribution of V± L \ S^ 2 J L) E n i, is nontrivial, and may be found in PYI. 
We write 



(3.3a) 



and the TV are given in PYI. These equations are unchanged from PYI; we repeat them here for ease of 
reference. 

The dominant nonperturbative corrections are associated with the gluon condensate and are' 4 ' 



x rp 2 / n 2\ f na \ A e n in 6 n{a s G 2 ) 

ONpE n i = me n m TT{a s G ) — = m— — - ; 

V 2 / (mOi?a s ) 4 

_ 1872 _ 2 102 _ 9 929 

10 1275' 20 1326' 21 9 945' 



(3.4) 



Because (a s G 2 ) ~ A 4 , this is of order (A/m) A . Besides the corrections reported above, there are a few 
pieces of the higher order perturbative and nonperturbative corrections that are known; they can be found 
discussed in PYI. 
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4 Numerical Results 

Using the formulas deduced above one evaluates the quark mass, and spectrum and other properties of 
heavy quarkonium systems. We take, 

A(n f = 4, three loops) = 0.23±g;g| GeV [ a s {M 2 z ) ~ 0.114±g;gg2 ] , (4.1a) 

and for the gluon condensate, very poorly known, 

(a s G 2 ) = 0.06 ± 0.02 GeV 4 . (4.1b) 

This value of a s (Af§) is slightly smaller than, though compatible with, the world average a s (M z ) = 0.118. 
We have preferred our value, which is obtained averaging measurements performed at spacelike momenta; 
see the recent review of S. Bethkel 8 !. 

Another matter to be discussed is the choice of the renormalization point, \i. As our equations (3.2, 3) 
show, a natural value for this parameter is 

li=—, (4.2) 
na 

for states with the principal quantum number n, and this will be our choice. For states with n = 1 the 
results of the calculation will turn out to depend little on the value of /j,, provided it is reasonably close to 
(4.2). Higher states are another matter; we will discuss our choices when we consider them. 

The 10 state of bb and the mass mj. As stated, we select, for the T state, /j, = 2/a. We then use Eqs. (3.1-4) 
to obtain the values of the b quark mass. The results are reported below; the errors correspond to the errors 
in Eqs. (4.1a, b). In the estimate of the errors, the condition \x = 2/a is maintained satisfied when varying 
A while for the error due to the variation of \i the other parameters are kept fixed (i.e., one no more has 
then n = 2/a). The dependence of mj on /j, should be taken as an indication of the theoretical uncertainty 
of our calculation. With Schroder's value for b\, 

m b = 5.001±g;g|I {A) =F 0.005 {(a s G 2 )) (vary ^ 2 by 25%) ± 0.006 (other th. uncertainty) 

to 6 (to 2 ) = 4.440±g;gi| ( A ) T °- 005 (( a sG 2 )) (vary ^ 2 by 25%) ± 0.005 (other th. uncertainty). ^'^ 

The values of yu 2 , a s (fi 2 ), a s (/j, 2 ) are, respectively, 

fj 2 = 6.632 GeV 2 , a^n 2 ) = 0.246 , a^u 2 ) = 0.386. 

The piece denoted by the expression "other th. uncertainty" in (4.3) refers to the error coming from 
higher dimensional operators and higher order perturbative terms; it can be found discussed in PYI. It is 
comfortably smaller that the errors due to the uncertainty on A, (a s G 2 ). If we omit these errors, so as not 
to double count them, and consider that the theoretical error is only that due to varying /j, 2 by 25%, and 
compose all the errors quadratically, then we obtain the estimates reported in the Introduction, Eqs. (1.1). 

M(Y) — M(rjb); the decay Y — > e + e~. The expressions for the hyperfine splitting, and the decay of the Y 
into e + e~ are as in PYI. They depend on b\ only indirectly, through the preferrred values of to, /i. We have 
the numerical results, using the Schroder calculation 

M(Y) - M( V ) = 46.6±J° 5 9 (A) j£« ((a s G 2 )) tf^ = 6-632 ± 25%) (4.4) 

and 

r(Y - e+e") = 1.07+8:oi -o.u (i^G 2 )) t° f 6 ( M 2 = 6.632 ± 25%). (4.5) 

practically unchanged from PYI. Note that, when varying A, {a s G 2 ) 7 we have varied nib according to 
Eq. (4.3), but we have not varied nib when varying \i. Note also that the corrections are here fairly large; in 
particular, due to the large size of the radiative correction to the decayJ 9 ' 

Higher order NP corrections due to the higher dimensional operators are also known for the decay rate 
(see PYI). Size corrections, however, are not known now. 
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The result for the decay is in reasonable agreement with experiment, 

r eX p.(T -> e+e") = 1.320 ±0.04 keV. 

Higher states (n = 2) of bb. As is clear from the expressions (3.2, 3) the natural choice of scale is now 
H = 1/a = 2.860 GeV 2 . If we take this, adding or subtracting a 25% to estimate the dependence of the 
calculation on the choice of scale then we obtain the results, with the Schroder value of b\ are, 

M(20, th) - M (20, exp) = 363^ MeV (fi 2 = 2.86 ± 25%), 

M(21, th) - M(21, exp) = 208±i?| MeV (/i 2 = 2.86 ± 25%). ' '"'' ! 

The 10 state of cc and the mass m c . The value of the parameter A used now, corresponding to that in 
Eq. (4.1a), is 

A(n f = 3, three loops) = 0.30±°;™ GeV . 
The values for the c quark mass, deduced from the J/ij) mass are now, 

m c = 1.8661^(^1) =p 0.014 ({a s G 2 }) (varying fi 2 by 25%) ± 0.011 (th. uncertainty) 

m c (m 2 ) = 1.531±g-g|| (A) =p 0.013 {{a s G 2 )) (varying fi 2 by 25%) ± 0.010 (th. uncertainty) (S), 

(4.7) 

and /i 2 = 2.465 GeV 2 now. Composing the errors as for the 6-quark we find the results reported in the 
Introduction, Eqs. (1.1). 

5 Discussion 

The discussion of PYI holds valid for the calculations using both Peter's and Schroder's evaluations, with 
one point of difference. If we believe Schroder's value of b\, then the two loop corrections are comfortably 
small. For example, with Peter's value we had 

a x a s /it ~ 0.11, M 2 /7r 2 ~ 0.14, (P) 

while with Schroder's the first is almost unchanged but the second becomes bia 2 /ir 2 ~ 0.081 (S) 
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